Hydrodynamics of a Micro-Hunter: Chemotactic Scenario 
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Inspired by bacterial chemotaxis we propose a hydrodynamic molecular scale hunter that can 
swim and find its target. The system is essentially a stochastic low Reynolds swimmer with ability 
to move in two dimensional space and sense the local value of the chemical concentration emitted by 
a target. We show that by adjusting the geometrical and dynamical variables of the swimmer we can 
always achieve a swimmer that can navigate and search for the region with higher concentration of 
a chemical emitted from a source. The system discussed here can also be considered as a theoretical 
framework for describing the bacterial chemotaxis. 

PACS numbers: 07.10.Cm, 87.17.Jj, 47.15. G- 



Propulsion mechanisms for microorganisms and artifi- 
cial swimmers are subject to the exceptional constraints 
of the motion in low Reynolds number hydrodynamics 
[H,Q. Purcell's Scallop theorem beautifully demonstrates 
how a set of non-reciprocal body deformations is neces- 
sary to achieve a net translational or rotational move- 
ments in simple systems (3l-[lo|. Experimental verifica- 
tion of the swimming motion in systems with only a small 
number of internal degrees of freedom have attracted in- 
terest in developping new artificial swimmers [ill [l2| . In 
biological systems, sperm cell as micron scale hunter, uses 
beating flagella to swim toward its targets These 
targets include food that is necessary for surviving and 
the egg cell that is essential in fertilization process. A 
concentration gradient of chemical emitted by the source 
is established in this chemotaxis phenomena [l7j]. The 
physical mechanism of chemotaxis in flagellated cells like 
sperm, with circling trajectories are usually described in 
the following way [18j : the underlying chcmichal network 
of an active stimulus-response system provides a concen- 
tration mediated stimulus which periodically regulates 
the internal motion and modulate the curvature of the 
swimming path. The input in this signaling system is 
the local value of a chemical 0, [2(| ■ This scenario gives 
rise to a drift in the circular trajectories of chiral flagel- 
lated swimmers 18l 21\. 

In this letter, inspired by chemotaxis we propose a 
molecular scale swimmer that can navigate in (2-D) 
space, along the gradient of a stimulating chemical. We 
assume that the chemical can activate a simple relax- 
ational process in the swimmer and correspondingly ini- 
tiates signals to change the internal motion leading the 
system to find the right track. The effects of fluctuations 
have been taken into account by considering a stochastic 
description of the internal motion of the swimmer. We in- 
vestigate the conditions that the swimmer can reach the 
region with higher concentration of the chemical. The 
model presented here can also be regarded as a theoreti- 
cal model to account the hydrodynamic details of chemo- 
taxis for a chiral swimmer with a single flagcllum. 

Consider a minimal hydrodynamic propeller composed 
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FIG. 1: (color online). Two large and small spheres are con- 
nected through an arm to construct a swimmer that can move 
in (2-D) space. This model resembles the geometry of a bac- 
terium that has a single tail. In a stochastic description of 
the motion we assume that the small sphere can be in one of 
the three distinct states as shown in the picture. The overall 
motion of the system in a complete cycle of the stochastic 
jumps are shown. 



of two large and small spheres with radii R and a where 
a < R. These spheres are connected through a negligible 
diameter arm that is not interacting with the ambient 
fluid. Different internal configurations of the system can 
be achieved by changing the length and the shape of the 
arm. As a specific choice we consider the case where the 
arm always stay in a two dimensional plane, as shown 
schematically in Fig. [TJ For further simplification of the 
model we decrease the number of internal configurations 
to three different states denoted by states (1), (2) and 
(3). State (1) is a reference state where two spheres are 
separated by a distance L. States (2) and (3) are char- 
acterized by a distance e and an angle $ with respect 
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to the reference state (1), as depicted in Fig. Q] The 
small sphere can jump between the denoted states along 
the shown paths with a constant velocity v Q measured 
in the reference frame of the large sphere. A complete 
cycle of the motion fulfills the Purcell's Scallop theorem 
and eventually leads the system to a new state with a 
net translational and rotational displacements Q. The 
motion is restricted to a (2-D) plane that is characterized 
by the plane of the arm. 

To analyze the hydrodynamics of the system, we use 
the fact that the linear Stokes equation governs the hy- 
drodynamics of the fluid flow at zero Reynolds number. 
To simplify the description, we furthere assume that the 
radius of small sphere is much smaller than the other 
length scales of the system, namely R and L (a R, 
a<L), In this case the motion of the small sphere can 
be regarded as a singular body force located at the posi- 
tion of this sphere. Taking the advantage of this simplifi- 
cation and using the linearity of the Stokes equation, we 
can describe the dynamics of the system for a general in- 
ternal motion. In the reference frame, that is co-moving 
and rotating with large sphere, we denote the position 
vector of the small sphere by ro(t). Denoting the linear 
and angular velocities of the large sphere in a laboratory 
frame by V and $7, we can express the velocity field of 
the fluid at a general point r in the co-moving frame as: 

u(r) = -V - n x r + M.V + ft x m + G(r, r ).f , (1) 



where f denotes the strength of the point force located 
at the position of the small sphere. Here the tensor M 
and vector m give the flow field due to the translational 
and rotational motion of a moving sphere and are given 
exphcitcly as: M = f f [1+ + 3£] and m = 

-p-r. The Green's function of the Stokes equation for an 
infinite flow bounded internally by a solid sphere with 
radius R is denoted by G(r, ro). The explicit form of this 
Green's function has been calculated by C.W. Oseen 22 1. 

Force and torque balances for a self propeller system 
require that the total force and torque acting on the fluid 
should vanish. The point force located near the sphere 
has an image force with strength f,; = (cifofo + C2I) • f 
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the unit matrix. This image force is located inside the 
sphere in position given by: r-Q = ^tTq. A discussion by 
Higdon shows that in evaluating the force acting on the 
fluid one should carefully account the image system 23 1 . 
In this case the force and torque balances read: 

f + fi + 6m]RV = 0, r x f + r* x f t + 87r?7i? 3 J7 = 0. 

To finish with the dynamical equations we should include 
the prescribed form of the internal motion by the bound- 
ary condition: u(ro) = io- Having in hand the force bal- 
ance equations and this boundary condition, we can elim- 
inate the point force strength and arrive at equations for 



linear and angular velocities of the system and obtain ve- 
locity of large sphere as: V = A • ro and O = B • ro , where 
A and B are two matrices with elements that strongly 
depend on the specific form of the internal motion given 
by function r (t) HJj ■ Here instead of giving the explicit 
expression of the swimming velocity for a general motion, 
we concentrate on the simple motion of our (2-D) model 
that has been introduced and discussed already. 

We denote the dynamical variables of the system by 
x and 9, where x stands for the position vector of the 
large sphere and 9 measure the angle that the swimmer's 
director makes with the x— axis. Swimmer's director is 
defined as a unit vector pointing from the position of the 
reference state (1) to the center of the large sphere. Note 
that x = V, 6 = H. Now the differential changes of the 
swimmer's variables in a general jump from state (i) to 
state (J) can be written as: 



Ax y - =H- 1 (9)-d lJ , A9 t 
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where R(#) represents the matrix for a rotation around 
z— axis by the instantanous value of the angle 9. Let 
us consider the case where the internal deformations are 
small, compared to the average length of the swimmers, 
such that e<i. This allows us to set up a perturbativc 
expansion of the results. Up to the leading order in e and 
a, the differential rotations read: 
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CJ23 = — 4ai2 and 031 = oti2- To the leading order in e 
and a, the displacement vectors read: 



di2 = 



Si 



123 



d 3 l = 



-81 
60 



where 



Si 



a , 
e-(l 



R 



3R 
2L 



)cos$, 62 



a , 



R 



3i? N 
IT' 



sin$, 



and 63 = — A82 and the results are given for R -C L. Scal- 
lop theorem allows us to simply express the changes for 
reverse jumps in terms of the forward jumps such that: 



Axji — —Axij and A£>,j = —A9ji. We have assumed 



that all jumps happen with a constant velocity vq. In 
this case the time for the jump 1 — > 2 (3 —> I) is equal to 
Ti = e/vo (r 3 = c/vq) and the time for the jump 2 — > 3 is 
equal to T2 = 2sin$e/i>o. The trajectory of the motion 
for the swimmer moving in a cyclic way (1 — >-2 — >-3 — >1) is a 
circular path with curvature given by: Kq 1 = (262 + S3). 
In this motion the swimmer's director rotates with an an- 
gular frequency given by: Oo = (2ai2 + ^23) / (2ri + r 2 ). 

To consider the effects of fluctuations which are present 
in any physical system, we develop a stochastic descrip- 
tion of the system, where the jumps between different 
states can happen with different rates. To construct the 



FIG. 2: (color online). Two different trajectories for the swim- 
mer moving in a linear gradient of concentration given by 
p(x) — a± + a,2X are shown. Parameters in both graphs are 
given by: R = 1, wo = 1, a = 0.1, L = 3.1, e = 0.4, a± — 
2, <J2 = 0.5, a = 1, /i = 10, for left graph: <1> = 0.2-7T, CJ13 = 
1.5 while for right graph: $ = 0.17T, ojis, = 1-5. For initial con- 
ditions we have chosen: a;(0) = 5, j/(0) = 0, 8(0) = 0. Two 
different behavior of the trajectories can be distinguished: 
motion toward higher or lower concentrations. 



the stochastic model we denote the probability of the 
system to be in the state (i) by Pi. The transition rate 
for stochastic jump from state (i) to state (j) is denoted 
by tUij. The rates for internal conformational changes 
in general depends on the temperature of the fluid and 
the detail internal activity of the system. Dynamics of 
this stochastic system is governed by the Fokkcr-Planck 
equations as: P\ = uj 31 P 3 + LO21P2 - {^12 +^is)Pi, Pi = 
<^32-P3 + w\iP\ — (c^2i + ^23)^2, where dot symbol de- 
notes the time derivative and also note that the proba- 
bility conservation implies that P\ + P2 + P3 = 1. The 
differential change of the displacement per unit time for 
this stochastic system can be written as: 
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Pi (W12AX12 + W13AX13) 



P 2 (W23AX23 + w 2 i Ax 2 i) + P 3 (w 3 i Ax 31 + W32AX32) , 

and a similar equation for the rate of change of 9 while 
we replace all Axy with oiij. For a system that is in 
the thermodynamic equilibrium, the transition rates for 
all jumps are symmetric (wjj = ujji) and the time aver- 
aged velocity of this system is zero. If for any reason the 
detailed ballance violates in the internal conformational 
changes, the system will behave like a circle swimmer. As 
an example let us assume that the rates for all the clock- 
wise jumps are equal, such that ui 12 = u 2 3 = w 3i = w o 
and all the other counterclockwise jumps set to co a + Scu. 
In this case and for small 5ui the trajectory is a circle with 
radius of curvature given by = (Suj/ujo)kq , where 
we have already defined kq. 

Now let us consider a fluid medium occupied by a very 
low concentration of chcmichal attractants given by p(x). 
The low concentration assumption makes sure that the 
hydrodynamic properties of the medium are not affected 
by this chemical. This chemical can derive the system to 
a nonequilibrium case by affecting the internal conforma- 
tional changes. Inspired by the chemotactic navigation 
in biological microorganisms, we assume that there is a 
density sensing mechanism in the system. In the sig- 
naling network of the chemotactic systems, a dynamical 



FIG. 3: (color online). Two different trajectories for the 
swimmer moving in a central gradient of concentration given 
by p(x) = 03/|x| are shown. Parameters for left graph are: 
$ = 0.2-7T, W13 = 0.8 and for right graph: $ = 2.6n, W13 = 
0.5, other parameters are the same parameters used in Fig. 
[2] Two different behavior of the trajectories can be distin- 
guished: motion toward higher or lower concentrations. 
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FIG. 4: (color online). Different behavior of the trajectories in 
a central concentration categorized in a phase diagram. The 
horizontal axis $ stands for a geometrical variable of the sys- 
tem and the vertical axis shows W13, the anisitropy parameter 
that makes the individual swimmer an active system. Spiral 
trajectories going toward (backward) the region with higher 
concentration where is located in the center, are denoted by 
A (R). Parameters have the same values as Fig. 2] 



mechanism is capable to produce an output signal that 
depends on the time history of the stimulating proper- 
ties. Local concentration of the chemotoattractants is 
the stimulating properties in this case. Mathematics of 
such a generic adaptation mechanism can be modeled by 
a simple relaxational process 0, [25[ . In this model an 
internal adaptation variable u(t) couples to a signaling 
function s(t) through the following equations: 



as 



pu 



pii = it(l — s), 



(4) 



where the adaptation variable u, measures the dynam- 
ical sensitivity with the time scale of adaptation given 
by p. The relaxational process in a time scale controlled 
by a will produce a stimulus s(t) that can affect the dy- 
namical variables of the conformational changes which 
are already defined by Uy. As an example we consider 
a case where the adaptation mechanism control only one 
of the transition rates. Here we choose W21 = s(£) w i2 
and set all the other rates including W12, to ujq. For a 
uniform profile of the concentration, the system reaches 
a steady state with s = 1 that is not moving. If we set 
an asymmetry in one of the other internal conformational 
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jumps, for example W13 ^ the steady state solution 
is a circling path. Now we can investigate the responcc 
of a circling swimmer, embedded in a gradient of chemi- 
cal concentration. First we consider a linear gradient of 
chemical concentration given by p(x) = a\ + C12X. Fig. 
[5] shows the trajectory of the swimmer for two different 
vales of the rates and the geometrical characteristics of 
the swimmer. The results can be summarized as follow: 
in general, nonuniform chemical concentration initiates 
a drift and the path looks like a drifting circle moving 
to the region with either higher or lower concentration 
value. For a constant asymmetry in one of the confor- 
mational changes, here W13 — W31 > 0, there is a critical 
value for $ c . Swimmer with geometrical structure that 
corresponds to <E> > <& c ($ < $ c ) move to the region with 
higher (lower) concentration. Interestingly all the results 
are independent of the initial orientation of the swimmer 
and are not sensitive to the adaptation paramctercs. 

Motion in central distribution of the concentration has 
the same characteristics as linear gradient. Fig. [3] shows 
the results for motion in a central gradient described by 
p(x) = a3/|x|, with 0:3 > 0. The trajectory in this case is 
a drifting circle along a spiral path. By changing the pa- 
rameters we can achieve a swimmer with spiral trajectory 
that is moving cither toward or away from the center. A 
phase diagram for different kinds of the possible states is 
shown in Fig. [5] In this phase diagram attractive states 
into the center (region with high concentration) and re- 
pulsive states from the center are denoted by labels (^4) 
and (R) respectively. The horizontal axis $, shows the 
geometrical variable of the system and the vertical axis 
(jJ±3, is the asymmetric parameter that makes the indi- 
vidual swimmer an active system. This picture shows 
that for an active swimmer (a system that can move in a 
uniform concentration), by adjusting the parameters it is 
always possible to find states which can navigate to the 
correct direction. 

In a concentration gradient, the force due to the 
nonequilibrium concentraton is also important [26| . Here 
we have neglected the effects due to this force. A di- 
mensional analysis shows that for satisfying this con- 
dition, the following criterion must hold: |Vp| <C 
u]o(aer))/(R 3 kBT), where luq is a typical value for the 
rates of conformational changes. In conclusion, we have 
introduced a micron scale system that uses the physics 
of chemotaxis in bacteria and navigate along a preferred 
direction into a source. The system is not sensitive to 
the adaptation variables. This two dimensional model 
can be readily generalized to three dimensional case. 
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